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Abstract 

We study the polytropic gas scenario as the unification of dark matter and dark energy. 
We fit the model parameters by using the latest observational data including type la su- 
pernovae, baryon acoustic oscillation, cosmic microwave background, and Hubble parameter 
data. At 68.3% and 95.4% confidence levels, we find the best fit values of the model param- 
eters as K = 0.742t^(iCT)tg;£g(2cr) and n = _i.05l°;°|(la)t°;^(2o-). Using the best fit 
values of the model, we obtain the evolutionary behaviors of the equation of state parameters 
of the polytropic gas model and dark energy, the deceleration parameter of the universe as 
well as the dimensionless density parameters of dark matter and dark energy. We conclude 
that in this model, the universe starts from the matter dominated epoch and approaches a 
de Sitter phase at late times, as expected. Also the universe begins to accelerate at redshift 
Zt = 0.74. Furthermore in contrary to the ACDM model, the cosmic coincidence problem is 
solved naturally in the polytropic gas scenario. 
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1 Introduction 



Various cosmological observations, including the type la supernovae (SNela) [1], cosmic mi- 
crowave background (CMB) [2] and baryon acoustic oscillation (BAO) [3], etc., have revealed 
that the current expansion of the universe is accelerating and it entered this accelerating phase 
only in the near past. A new energy component with negative pressure called dark energy (DE) 
is needed to explain this acceleration expansion within the framework of general relativity (GR). 
The simplest and most appealing candidate for DE is the A cold dark matter (ACDM) model, in 
which the cosmological constant A plays a role of DE in GR. This model is in general agreement 
with current astronomical observations, but has difficulties in reconciling the small observational 
value of DE density to that coming from quantum field theories. This is called the cosmological 
constant problem [4]. Therefore, a number of models for DE to explain the current acceleration 
without the cosmological constant has been proposed, such as quintessence [5], phantom [6], 
K-essence [7], tachyon [8], quintom [9]; as well as the Chaplygin gas [10] and the generalized 
Chaplygin gas [11], the holographic DE [12], the new agegraphic DE [13], the Ricci DE [14], and 
so on. 

One of interesting DE models is the polytropic gas which was proposed to explain the accel- 
erated expansion of the universe [15, 16, 17]. It was shown that the polytropic gas model in the 
presence of interaction with DM can behave as phantom type DE [16]. It was pointed out that a 
polytropic scalar field can be reconstructed according to the evolutionary behaviors of the holo- 
graphic [18] and new agegraphic [19] DE densities. The validity of the generalized second law 
of gravitational thermodynamics was examined for the polytropic gas model in [20]. Different 
modified gravity models including /(T)-gravity [21] and /(i?)-gravity [22] were reconstructed 
according to the polytropic gas equation of state. The polytropic tachyon, K-essence and dilaton 
scalar field models were studied in [23] . In [24] , the polytropic gas model was investigated from 
the viewpoint of statefinder diagnostic tool and w — w' analysis. 

In the present work, our main aim is to fit the polytropic gas model and give the constraints 
on model parameters, with current observational data including SNela, CMB, BAO and obser- 
vational Hubble data (OHD), filling in the gap existing in the literature. To do so, in section 
2, we briefly review the polytropic gas model as unified picture of DM and DE in a spatially 
flat Friedmann-Robertson- Walker (FRW) universe. In section 3, we consider the cosmological 
constraints on the polytropic gas model by using the latest observational data. In section 4, we 
give numerical results. Section 5 is devoted to conclusions. 

2 Polytropic gas model 



For a polytropic gas model, the energy density p po \ and pressure p po \ satisfy the equation of 
state (EoS) [17] 



where K > and the polytropic index n are two constants of the model. 

Here, we consider a spatially flat FRW universe containing only the polytropic fluid and 
baryonic matter. We ignore the contribution of radiation, as expected from the observations 
[25]. In the framework of the standard FRW cosmology, the first Friedmann equation reads 



where H = a/a is the Hubble parameter. Also p po \ and /3bm denote the energy densities of the 
polytropic and pressureless baryonic matter (pbm = 0), respectively, which satisfy the following 



P P o\ = - K P P oi 



(1) 




(Ppo\ + />bm), 



(2) 
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energy conservation laws 

Ppol + 3H(p po i + pp i) = 0, 
Pbm + 3Hp hm = 0. 
Using the dimensionless density parameters 

_ Ppol _ 87rGp po i _ Pbm 

Per oH * p cr 

the Friedmann equation (2) is rewritten as follows 

1 — ^pol ~i~ ^bm' 



8-7rGpbm 

3H 2 ' 



(3) 
(4) 

(5) 
(6) 



Substituting Eq. (1) into (2) gives the evolution of the polytropic gas energy density in terms 

(7) 



of redshift z = - — 1 as 

a 



Ppol Ppolo 



where 



K+ (1 -K)(l + z)- 3/n 

l/n 



K = P P ol 



K, 



and p P oi is the polytropic energy density at present time. Also Eq. (4) gives 

Pbm = Pbm (l + zf. 



(8) 
(9) 



It is worth to note that the polytropic gas energy density (7) offers a unified picture of DM and 
DE. Because it smoothly interpolates between a non-relativistic matter phase, p po \ oc (1 + z) 3 , 
in the past and a negative-pressure DE regime, p po \ = —p po \, at late time (z — > —1). Therefore, 
within the framework of unified DM-DE (UDME) scenario one can rewrite the polytropic gas 
energy density as 



Ppol — Pdm + Pdo 
Ppol = Pdc, 



(10) 



where the DM is a pressureless matter (pdm = 0). Using Eqs. (1) and (7) one can obtain the 
EoS parameter a; po i of the polytropic gas model as 



Wpol(z) 



Ppol ^ 
Ppol 



K 



K + {1-K){1 



-3/n ' 



(11) 



With the help of Eqs. (3), (4), (7) and (10), the energy densities of DM and DE evolve as follows 



Pdm = Pdm (l +zf, 



Pdc = Ppolo 



K + (l-K)(l + z) 



-3/n 



- Pdm (l + zf. 



Using Eqs. (1), (5) and (13) the EoS parameter of DE reads 

Pde 



Ude{z) = 



Pdc 



-k(i 


^bmo ) 


k + (1- K)(l + z)-^/ n 


-1-n 


(1 - fibmo) 


k + (i- 


- k)(i + z)- 3 / n 


—n 

-f] dmo (l + z) 3 



(12) 
(13) 

(14) 
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With the help of Eqs. (7) and (9) and using (5), the first Friedmann equation (2) takes the form 



H 2 (z,p) = H^(l - n hmo )[K + (I - K)(l + z)- 3/n ] n + fWl + ^) 3 }, (15) 



where H = 70.2±1.4 km s" 1 Mpc" 1 (68% CL) and Q hmo 



8nGp hn 
3H% 



0.0458±0.0016 (68% CL) 



are the present Hubble constant and the present value of the dimensionless BM density, respec- 
tively, which have been updated in the 7-year WMAP (WMAP7) data [25]. Also p indicate 
model parameters including K and n. Thus, throughout this work we fix the Hubble and the 
baryon density parameters at Hq = 70.2 and f2bm = 0.0458. With H$ and rit> mo being deter- 
mined by independent measurements, in the next section we will use the cosmic observations to 
constrain the polytropic gas model parameters (K, n). 
For completeness, we give the deceleration parameter 



q=-l- 



H 



(16) 



which combined with the EoS and the dimensionless density parameters form a set of useful 
parameters for the description of the astrophysical observations. Using Eq. (15) the deceleration 
parameter (16) can be obtained as 



^bm (l + Z? + (3^pol + 1)(1 - ^bm ) 


k + (i- k)(i + z)- 3 / n 


—n 


^bm (l + ^) 3 + (l-^bm ) 


k + (i- k)(i + «)-3/« 


—n 



(17) 



3 Observational constraints 



Here, we fit the free parameters of the polytropic gas model (1) by using the recent observational 
data including SNela, BAO, CMB and OHD. 

For the SNela data, we use the currently largest Union2.1 compilation [26] that contains a 
total of 580 SNela, which is an updated version of the Union2 compilation [27]. Cosmological 
constraints from SNela data are obtained through the distance modulus n(z). The theoretical 
distance modulus is defined as [28, 29] 



fhh(z) = 51og 10 -D L (z) + 



(18) 



where hq = 42.38 — 51og 10 /i and h is the Hubble constant Hq in units of 100 km s 1 Mpc 1 . 
Also the Hubble- free luminosity distance Dl{z) for the flat universe is given by 



d L (z) = (i+z) r 

Jo 



dz' 



E(z';pY 



(19) 



with p the model parameters and E(z; p) = H(z; p)/Hq. 

For using SNela data, theoretical model parameters are determined by minimizing the quan- 
tity [28, 29] 

Xsn = A-^, (20) 



where 



580 



(21) 



i=l 
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580 

B = £[/W*0 " fM^(zi)]/of, (22) 

i=l 

580 

^ = E!M 2 , (23) 

i=l 

and <7i stands for the la uncertainty associated to the ith data point. 

Next, we add the data from the observation of acoustic signatures in the large scale clustering 
of galaxies. Using the BAO data, one can minimize the Xbao defined as [30, 3], 

Xbao = KbS ~/ th]2 , (24) 

where 



A th = JnZ E^py 1 / 3 - f 

V l z b JO 



1 f z b dz' 



2/3 

(25) 



_z h J E(z';p)_ 

is the theoretical distance parameter and z^ = 0.35 is the redshift of luminous red galaxies 
sample of the Sloan Digital Sky Survey (SDSS). Here A ohs = 0.469(n s /0.98)~ a35 ± 0.017 is 
measured from the SDSS data [3] and the scalar spectral index n s is taken to be 0.968 [25]. 

Since the SNela and BAO data contain information about the universe at relatively low 
redshifts, we will include the CMB shift information by using the WMAP7 data [25] to probe 
the entire expansion history up to the last scattering surface. The shift parameter R of the 
CMB is defined as [31, 32] 

where z rec ~ 1091.3 is the redshift at the recombination epoch [25]. Also O mo is the effective 
matter density parameter defined as 

ftmo = ^bmo + (1 " flbmoXl " K)~ n . (27) 

This expression for Q mo , is an estimate of the "matter" component of the polytropic gas fluid 
with the baryon density. The x 2 from the CMB constraint is given by 

2 _ ffiobs - -Rth] 2 

where the observational value of i? Q b s has been updated to 1.725 ± 0.018 from the WMAP7 data 
[25]. 

Finally, we further add the data from the observational Hubble parameter. The Hubble 
parameter is related to redshift with 

HU) = (29, 

so if dz/dt is known, H(z) is obtained directly [33]. The best fit values of the model parameters 
from OHD can be determined by minimizing [34] 

„2 _ ST [Hobs(Zi) - Ht^Zj^)] 2 

XOHD - 2^ Z2 ' 

i=l i 
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where is the theoretical value of the Hubble parameter and H ^ s is the observed value. Table 
1 shows the observational Hubble data containing the nine data from [35] and three additional 
data, in bold face, from [36]. 

As the relative likelihood function is defined by C = e~ ( - x totai _ x m i n )/2 ^S7], the best fit value 
of the model parameters follows from minimizing the sum 

Xtotal = XSN + XBAO + XCMB + XOHD- (31) 

4 Numerical results 

The best fit values and errors of the polytropic gas model parameters (K, n) are summa- 
rized in Table 2, where we also list the best fit results of the ACDM model for compari- 
son. At la (68.3%) and 2a (95.4%) confidence levels (CLs), we obtain the best fit values 
K = 0.742^1| (la)tox>il(2<r) and n = — 1 -05±g;g| ( l<r)±g;i| (2a) for the full data sets includ- 
ing SNela+BAO+CMB+OHD. The total \ 2 of the best fit values of the polytropic gas model 
is Xmin = 573.089 for the full data sets with dof (degree of freedom) = 594. The reduced 
X 2 = Xmin/ 'd°f 1S 0-965, which is acceptable. The obtained Xmin ^ s slightly smaller than the 
one for the ACDM model, Xacdm = 573.552, for the same data sets. The marginalized relative 
likelihood functions C{K) and C{n) are shown in Figs. 1 and 2, respectively. Considering only 
the SNela data, the degeneracy between the model parameters is considerable, but using the 
combined data this degeneracy is dropped in different ranges of redshift. Figure 3 shows the 
constraint on the polytropic gas parameter space n — K at la and 2a CLs, using the full data 
sets. 

The variations of the EoS parameters of the polytropic gas model u; po i and DE Wdc with the 
best fit values of the model are plotted in Figs. 4 and 5, respectively. Figure 4 presents that 
Wpoi at early (z » 1) and late (z — > — 1) times behave like the EoS parameters of the matter 
(wpoi = 0) and the cosmological constant (oj po \ = —1), respectively. This shows that in the 
polytropic gas scenario, the DM and DE can be unified. Figure 5 shows that Wde varies from 
Wde > — 1 to = — 1, which is similar to the freezing quintessence model [38]. The current best 
fit value of the EoS parameter of the DE in the polytropic gas model is obtained as Wdeo = —0.98 
which is in good agreement with the recent observational result Ude = —0.93 ± 0.13 (68% CL) 
deduced from the WMAP7 data [25]. 

In Fig. 6, we plot the evolutionary behavior of the deceleration parameter of the universe 
with the best fit values of the polytropic gas model, Eq. (17), and the ACDM model. Figure 
6 shows that very similar to the ACDM model the universe transits from an early matter 
dominant regime, i.e. q = 0.5, to the de Sitter phase, i.e. q = —1, in the future, as expected. 
The accelerating expansion begins at transition redshift z t = 0.74, which is earlier than what the 
ACDM model predicts, ^ CDM = 0.72. The current best fit value of the deceleration parameter 
in the polytropic gas model is obtained as qo = —0.56 which indicates the expansion rhythm 
of the current universe. This is in good agreement with the recent observational constraint 
q = -0.43i£;}f (68% CL) obtained by the cosmography [39]. 

In Fig. 7, we plot the evolutionary behaviors of the energy density parameters of DM, 
!] dm = 87r 3 C ^ dm , and DE, £l de = 8 g^ dc , with the best fit values of the polytropic gas model 
using the full data sets. Figure 7 shows that fidm an d f^e decreases and increases, respectively, 
during history of the universe. We also obtain f2dm = 0.229 and f2d eo = 0.725 as the current 
best fit values. These are in exact agreement with the latest observational results f2dm = 
0.229 ± 0.015 (68% CL) and fl dco = 0.725 ± 0.016 (68% CL) deduced from the WMAP7 data 
[25]. We also get ^deo/^dmo = 3.166 ~ 0(1). This shows that in contrary to the ACDM model, 
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the cosmic coincidence problem, namely why the ratio of the DE and DM densities is of order 
unity today, is solved naturally in the polytropic gas scenario. 

5 Conclusions 

Using the latest observational data from SNela, BAO, CMB and OHD, we fitted the param- 
eters of the polytropic gas model as the unification of DM and DE. We obtained the con- 
straint results of polytropic gas model parameters, K = 0.7421q'q24(1c)^o'o49(2o') and n = 
-1.05to;o|(lcr)t£;^(2<7) for the full data sets. The minimal x 2 gives xLm = 573.089 with 
dof = 594. The reduced x 2 equals to 0.965 which is acceptable. The Xmm * s slightly smaller 
than the one for the ACDM model, Xacdm = 573.552, for the same data sets. 

Using the best fit values of the polytropic gas model parameters, we also studied the evo- 
lutionary behaviors of the EoS parameters of polytropic gas model and DE, the deceleration 
parameter of the universe as well as the DM and DE dimensionless density parameters. Our 
numerical results show the following. 

(i) The evolutionary behavior of the EoS parameter oj po \ of the polytropic gas model shows 
that the universe transits from an early matter dominated phase, i.e. oj po i = 0, to the ACDM 
model, i.e. u po i = —1, in the future, as expected. This confirms that the polytropic gas model 
plays the role of a unified model for DM and DE. 

(ii) The EoS parameter uj^c of DE in the polytropic gas scenario varies from w^c > — 1 to 
Wde = — 1 like a freezing quintessence model. The present value of w<j eo = —0.98 is in good 
agreement with the result of WMAP7. 

(iii) The variation of the deceleration parameter q shows that the universe transits from an 
early matter dominant epoch, i.e. q = 0.5, to the de Sitter era, i.e. q = — 1, in the future, 
as expected. The accelerating expansion begins at Zt = 0.74, which is earlier than that of the 
ACDM model. The current value of qo = —0.56 is in good accordance with the constraint coming 
from the cosmography. 

(iv) The evolutions of the dimensionless density parameters of DM and DE clear that f^dm 
decreases and f^e increases during history of the universe. The present values of dimensionless 
DM and DE densities are obtained as f2dm = 0.229 and £lde = 0.725 which are in exact 
agreement, surprisingly, with the results of WMAP7. Also the obtained present density ratio 
fideo/^dmo = 3.166 ~ 0(1) shows that in contrary to the ACDM model, the cosmic coincidence 
problem is solved naturally in the polytropic gas scenario. 
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Table 1: The observational H(z) data [35, 36]. 



z 0.09 0.17 0.24 0.27 0.34 0.40 0.43 0.88 1.30 1.43 1.53 1.75 



H(z) 


69 


83 


79.69 


70 


83.80 


87 


86.45 


117 


168 


177 


140 


202 


la 


±12 


±8.3 


±2.32 


±14 


±2.96 


±17.4 


±3.27 


±23.4 


±13.4 


±14.2 


±14 


±40.4 



Table 2: The best fit values of the parameters K and n within the 68.3% (la) and 95.4% (2a) 
CLs for each observational data set for the polytropic gas model. Columns 5, 6, and 7 show 
the current EoS parameter of DE, the current deceleration parameter of the universe and the 
transition redshift, respectively. The last row shows the best fit result of the ACDM model using 
the full data sets for comparison. 



Data 


K 


n 


x 2 - 

A mm 


^de 


<7o 




SN 


n 7c;o+0.031+0.055 
U./OO_ .034-0 .063 


-1.04 


562.225 


-0.98 


-0.59 


0.77 


SN + BAO 


n 7«n+°- 030 + - 057 
u -' DU -0.031-0.063 


i m +0.11+0.21 
— i - ui -0.12-0.25 


562.229 


-1 


-0.59 


0.75 


SN±BAO 
+CMB 


D 7^/1+0.028+0.054 
U./O^_ .028-0 .057 


i n o+0.08+0.15 

J^-^-o.og-o.is 


562.366 


-0.98 


-0.58 


0.75 


SN±BAO+ 
CMB±OHD 


fl 7/10+0.024+0.048 
u - '^-0.024-0.049 


i n c:+0.08+0.15 

J^-^-o.os-o.ie 


573.089 


-0.98 


-0.56 


0.74 


ACDM 






573.552 


-1 


-0.58 


0.72 
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Polytropic gas model 





Figure 3: The 68.3% (la) and 95.4% (2a) CL contours for n versus K from the full data sets. 
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Polytropic gas model 




-1 1 2 3 4 5 



Figure 4: The best fit of the EoS parameter uj po \ of the polytropic gas model, Eq. (11), using 
the full data sets. 
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Figure 6: The best fit of the deceleration parameter q{z) of the universe for the polytropic gas 
model, Eq. (17), and the ACDM model using the full data sets. 
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